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Abstract. Weak lensing surveys provide the possibility of identifying dark matter halos based on their total matter content 
rather than just the luminous matter content. On the basis of two sets of observations carried out with the CFHT, Erben et al. 
(2000) presented the first candidate dark clump, i.e. a dark matter concentration identified by its significant weak lensing signal 
without a corresponding galaxy overdensity or X-ray emission. 

We present a set of HST mosaic observations which confirms the presence of an alignment signal at the dark clump position. 
The signal strength, however, is weaker than in the ground-based data. It is therefore still unclear whether the signal is caused 
by a lensing mass or is just a chance alignment. We also present Chandra observations of the dark clump, which fail to reveal 
any significant extended emission. 

A comparison of the ellipticity measurements from the space-based HST data and the ground-based CFHT data shows a 
remarkable agreement on average, demonstrating that weak lensing studies from high-quality ground-based observations yield 
reliable results. 
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1. Introduction 

In the currently favored cosmological model, structure forma- 
tion in the universe is dominated by collisionless Cold Dark 
Matter (CDM). The model of structure formation by gravita- 
tional collapse in a pressure-less fluid is able to successfully 
reproduce the filamentary lar ge scale structure observed in the 
universe (e.g. |Peacocklll999l) However, for the formation of 
galaxies and galaxy clusters, gas dynamics play an important 
role. It seems obvious that galaxy formation is triggered when 
gas falls into the potential wells of dark matter concentrations. 
We therefore expect to find galaxies at the high-density peaks 
of the dark matter distribution. In the CDM scenario, small ha- 
los collapse earlier and merge to larger halos later. For galaxy 
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formation, this implies that the large galaxies we see today 
formed from mergers of protogalaxies. Observations support 
this theory: we see more irregular, small galaxies at higher 
redshifts and many merger systems and galaxies showing ev- 
idence for recent mergers. This bottom-up scenario also calls 
for galaxy clusters to build up through the merger of smaller 
halos. 

While it may be possible to (temporarily) drive gas from 
galaxy-sized halos, when such halos merge to form cluster- 
sized objects, the majority of them should contain galaxies 
and/or hot gas, so that the resulting massive halo is expected 
to contain a substantial amount of luminous matter. A cluster- 
sized halo very poor of luminous matter {dark clump) would 
require a mechanism to drive the gas out of all the smaller ha- 
los from which it assembled or from the massive halo itself. 
Both cases are highly unlikely: the first is very improbable, the 
second very difficult due to the high mass of the object. At the 
moment, there are no well-motivated physical processes to ex- 
plain either scenario. 

The discovery of a dark clump would therefore call for a 
critical reevaluation of our current understanding of structure 
formation in the universe. Currently, the only tool available to 
search for such objects is gravitational lensing, as it probes mat- 
ter concentrations independent of their nature. 

In the course of a weak lensing study. lErben et alJ |2000) 
announced the possible discovery of a dark clump, about 7' 
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south of the galaxy cluster Abell 1942. This assertion is based 
on significant alignment signals seen in two independent high- 
quality images, taken with the MOCAM and UH8K cameras at 
the CFHT. There is no associated apparent galaxy overdensity 
visible in these im ages nor in deep H-band images analyzed by 
iGrav et all fcOOll) . There is faint X -ray emission about I' from 
the lensing centroid detected by the ROSAT survey, but it is 
unclear whether this could be associated with a lensing object. 
If the alignment signal is due to a lensing mass at a similar 
redshift as the cluster, z = 0.223, this halo would have a mass 
of the order of 1O 14 /j _1 M . At a higher redshift (0.8 - 1.0), it 
would require a mass of the order of lO l5 h^ l M e . 

There are currently three more such dark clump candidates 
in the literature: 

- lUmetsu & Futamasel |2P00) find a candidate in their weak 
lensing analysis of the galaxy cluster CL 1604+4304 using 
data from the WFPC2 camera of the HST. In two separate 
datasets, they find a peak I'.l southwest from the cluster 
center, which corresponds to about 830 h~ l kpc at the red- 
shift of the cluster (z = 0.897). They estimate the mass of 
the object to be about 4.8 x 1O 14 /i _1 M , assuming it is lo- 
cated at a similar red shift as the cluster. 

- Miralles et al 1 ll2002l) found a conspicuous tangential align- 
ment of galaxies in an image taken by the STIS camera 
aboard the HST as a parallel observation. However, follow- 
up wide-field observations with the VLT failed to detect a 
weak lensing signal (lErben et alJl2003l) . so that a chance 
alignment of 52 galaxies in the original STIS analysis is at 
this point considered the most plausible explanation for this 
candidate. 

- IPahleetaTl 12003) identify a dark clump candidate about 
6' southwest of the galaxy cluster Abell 959 (z = 0.286) 
in images taken with the UH8K camera at the CFHT with 
evidence that this is a dark sub-clump of the cluster. If this 
is indeed an object at the redshift of the cluster, they deduce 
a mass of (1.1 + 0.3) x 10 1 4 h~ l M Q . 

IWeinberg & Kamionkowskil (2002) argue that about one 
out of five Dark Matter halos identified by weak lensing should 
be a non-virialized halo, i.e. a halo which is in the process 
of collapsing and has not yet reached dynamical equilibrium. 
Such objects should have only very little X-ray emission and 
about half the projected galaxy density of virialized halos. The 
luminosity of such objects would therefore be very difficult 
to determine. Accordingly, distinguishing between pure Dark 
Matter halos and normal, non-virialized halos may be almost 
impossible in these cases. 

However, the noise in weak lensing analyses due to intrin- 
sic ellipticities of galaxies can have a profound effect on the 
statistics of the number of halos detected per area. Intrinsic 
ellipticities may mimic tangential alignment, thereby causing 
fals e peaks or boosting the significance of lensing signals (e.g. 
iHamanaet all20ol. 

To determine the nature of the dark clump near Abell 1942, 
we obtained a set of HST observations of the field (General 
Observer Program, Proposal ID 9132, PI Erben). The HST 
probes fainter and thus more distant galaxies, for which the 



distortions of a foreground lensing mass are larger; addition- 
ally, due to the lack of seeing, its shape measurements should 
be more reliable. If the alignment signal seen in the ground- 
based data is due to a lensing mass, it should thus be even more 
significant in the HST data. 

The structure of this paper is as follows. In Sect. |2] we 
present the optical data available to us, namely the CFHT im- 
ages of the original detection and the HST data, along with 
our data reduction methods to extract object catalogs suitable 
for lensing studies. Sect.|3]gives a brief overview of the weak 
lensing methods employed in this paper. In Sect. 0]we present 
a re-analysis of the /-band image of the CFHT data. Our weak 
lensing analysis of the HST data, which confirms the alignment 
signal, but not its strength, is described in Sect. [5] In Sect.|6]we 
use a deep Chandra image to show that the ROSAT source is 
likely to be a spurious detection. The appendices illustrate var- 
ious tests for systematics (Adp.IaI and a comparison of the in- 
dividual shape measurements of objects common to the CFHT 
and HST datasets (App.|S}. 

2. Optical data 

The goal of th is work is to understand the origin of the lensing 
signal seen bv lErben etall koOO). We therefore consider both 
the ground-based dataset used in the original discovery as well 
as the HST data. Such a treatment also allows for a direct com- 
parison of the ellipticity measurements of objects detected in 
both datasets. 

2.1. Ground-based data 

Our anal ysis concen trates on the same /-band image as used 
in lErben e t al. (2000), as this covers most of the area imaged 
by the HST. We use Chip 3 of a mosaic observation of Abell 
1942 taken with the UH8K camera, with a pixel scale of 0'.'206. 
9 exposures of 1200 s went into the final image, which has 
a seeing of C'74. Unfortunately, a photometric calibration is 
missing. 

2.2. Space-based data 

Our HST data is a WFPC2 mosaic (approximately 5' x 4') of 
six pointings, each consisting of 12 dithered exposures with an 
exposure time of 400s each, taken between May 20th and June 
1st, 2001. The position of the mosaic with respect to the /-band 
image from the CFHT is shown in Fig.[^The filter employed 
wasF702W. 

2.2.1 . Data reduction 

Our reduction of the HST d ata is based largely on the dither 
package (tFruchter & Hookl l2002l) for IRAF. Each of the six 
pointings was reduced separately. Simultaneous processing of 
the four individual chips is done automatically by the dither 
routines. Due to their better signal-to-noise behavior, only the 
chips of the Wide Field Camera, namely Chips 2, 3, and 4 of 
WFPC2, were used for the later analysis. The dither pattern of 
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Fig. 1. The outlines of the six HST pointings superposed on the 
CFHT /-band image. The center of Abell 1942 lies in the top 
half of the CFHT image, the HST mosaic is centered on the 
dark clump position. 



the images allows us to achieve a higher resolution in the coad- 
ded image via the drizzle algorithm jGonzaea et all l998). 
The steps of the reduction are outlined in the following: 

Rough cosmic ray removal: To find the offsets between the im- 
ages, they first have to be cleaned of cosmic rays, which 
would otherwise falsify a cross-correlation. Each frame is 
cleaned using the precor task. This leaves only objects of a 
minimal size in the image, which should be stars and galax- 
ies, with little contamination by cosmic rays. 

Offset estimation: As the individual frames are dithered with 
respect to each other, it is necessary to find their relative 
offsets. This is done by performing a cross correlation of 
the cosmic ray cleaned images produced in the previous 
step. 



Median coaddition: Using the previously determined offsets, 
the original images are median combined. They are mapped 
via the drizzle algorithm onto an output grid with pixels 
of half the original size. 

Mask creation: The median image is mapped back onto the 
original resolution and offset of the original frames using 
the blot algorithm, the inverse of drizzle. The original 
frame is then compared to the median image to identify cos- 
mic rays via the deriv and driz_cr tasks. Thus, for each 
frame a mask is created identifying cosmic rays. This is 
combined with a mask identifying defect or possibly prob- 
lematic pixels, which is supplied with each raw frame. 

Offset determination: For each frame, those pixels that are 
flagged in the mask are substituted by their value in the 
blotted image (the transformed median image). These im- 
ages are then cross-correlated to determine the offsets more 
precisely. Possible small rotation angles have to be found 
manually. 

Coaddition: The images are drizzled onto an output grid of 
half the pixel size (pix-frac = 0.5), i.e. twice the origi- 
nal resolution. The value of each output pixel is obtained 
via averaging, where pixels which are flagged in the mask 
image are omitted. The drop-size (i.e. a scaling applie to 
the input pixels before being mapped) used is 0.6 . 

Performing this routine on all of the six pointings, we obtain 
18 single-chip images. A mosaic of these is shown in Figure|2] 

2.3. Catalog extraction 

For both datasets, we used a similar process to build the object 
catalogs. Differences in the procedure arise mainly from the 
small field-of-view and mosaic nature of the HST data. 

2.3.1 . Preparatory steps 

We manually updated the astrometric information of the 
ground-based image until its bright objects coincide with thei r 
positions as given in the USNO-A2 catalog jMonet etalJl 998). 
This step is not necessary for the lensing analysis of the ground- 
based image, but provides us with a reference catalog for the 
astrometric calibration of the HST mosaic. 

The HST images are aligned roughly with sky coordinates, 
but since each chip is read out along a different chip border, 
the individual images are rotated with respect to each other. To 
avoid confusion, we first rotate the images of Chips 2 and 4 
by +90°, respectively, such that north is up and east to the left 
(approximately). This is performed solely as a rearrangement 
of pixel values to avoid any additional resampling process. 

2.3.2. Masking 

In order to avoid noise signals and distorted ellipticity mea- 
surements, we mask out bright stars and the artefacts they cause 
(diffraction spikes, blooming, CTE trails, a ripple-like structure 
at the eastern edge of the /-band image, and brightening of the 
background along columns of the HST chips in two cases). In 
the HST image, we also mask the bright galaxy in the field. 
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Fig. 2. The HST mosaic with an illustration of the v arious dark c l ump c entroids cited in the text: the centroid we found in the 
ground-based data ( * ); the original position given bv lErbenetaDfeOOOl) ( + ); the peak position in the HST data ( A ); the center 
of the X-ray emission as found by ROSAT ( * ) and by Chandra ( x ); the position of the peak found in the 200" filter scale 
for the faint ground-based galaxies ( o ); and the galaxy number overdensity of medium-bright HST galaxies ( O ). Additionally, 
we plot a circle of 120" radius around the first position. For our M ap filter function of this radius, the maximal weight is then 
assigned along the dashed circle. Assuming a typical uncertainty of 1' for centroids found via weak lensing, all cited position are 
compatible with the dark clump. 



2.3.3. Source extraction 

We used SExtractor (Bert in & Arnoutsl 19961) to identify ob- 
jects in the image. SExtractor considers N connected pixels 
that are at a level k- <x sky above the sky background as an object, 
with cr sky being the standard deviation of the background noise. 
For the CFHT image, we used N = 3 and k = 1.0. These are 
very low thresholds, but since we later want to correlate objects 
present in both datasets, we strive to obtain a high number den- 
sity of objects. For the HST images we used N - 3 and k = 1 .5. 

HST: astrometric and photometric calibration 



The astrometric calibration plays an important role for a 
mosaic dataset such as our HST images, as it gives the position 
of the images with respect to each other. Indeed, the positions 
on the sky will be used later on in the lensing analysis rather 
than x and y position on the chip. 

For a reference catalog, we had extracted a catalog of bright 
objects from the astrometrically calibrated ground-based I- 
band image (see Sect. l2.3.T1 >. We match the entries of the refer- 
ence catalog to objects found in the HST images. This allows 
the determination of the pointing and the distortion of the im- 
age. Objects that are detected both in the /-band and the HST 
image can then be identified by sky coordinates. 
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The photometric calibration is done based on the relevant 
keywords of the HST image headers. 

2.3.4. Ellipticity measurement 

For the objects in the SExtracted catalogs, we measure the el- 
lipticities using a modifi ed version of the i mcat software, fol- 
lowing the method of iKaiser et alJ ([1995). We used the half- 
light radius as measured by SExtractor as the radius of the 
weight function with which the brightness profile is weighted. 

The measured image ellipticity \ is related to the source 
ellipticity x {s) by 



x =x (s) + P g g + P sm q* 



(1) 



where g is the reduced shear induced by a lensing mass which 
we ultimately want to determine, q* is the stellar anisotropy 
kernel, i.e. the anisotropy of the PSF for point-like objects. 
P g and P sm (smear polarizibility tensor) are tensors describ- 
ing a galaxy's susceptibility to the two distorting effects. 
The y are also measured from a galaxy's light distribution (see 
Bartelma nn & Schn eider 2001, for a derivation). 



2.3.5. Anisotropy correction 

To correct for the anisotropy induced by the telescope-detector 
system, we measure the ellipticities of the stars present in the 
field; for these eq. Q simplifies to 



(2) 



Stars are selected from a magnitude vs. radius plot. We fit a 
third-order polynomial in chip position (x, y) to the quantity 

( 0.5 *p»") to estimate Q*<X.y)- 

Thus, we obtain anisotropy-corrected ellipticity measure- 
ments: 



: X - P SI V 



(3) 



2.3.6. Modification for the HST images 



Since the PSF cannot be described analytically across chip- 
borders, the anisotropy correction for a mosaic has to be ap- 
plied to the single chips. With the small field of view of the 
WF chips, we have the added difficulty that for each image, 
there are only about five stars that could be used for the polyno- 
mial fitting, obviously not enough. Since the images were taken 
consecutively, we can assume that the PSF does not change 
considerably between the six pointings. We therefore apply an 
anisotropy correction for each chip based upon all the stars that 
were found in the six images taken by that chip. 

The 18 single-chip catalogs are combined to three catalogs, 
one for each chip. Because there are still few stars even in these 
catalogs, the stellar sequence is selected manually. For each 
catalog, a third-order polynomial is fitted. 

The anisotropy correction of the HST images is further dis- 
cussed in Arjpendix lA.il 



2.3.7. Shear estimation 

After the anisotropy correction, the second step in retrieving an 
estimate of the local shear from ellipticity measurements is the 
correction for the P s tensor. It is a combination of P sm and the 
shear polarizibility tensor P sh : 

pg _ psh _ psm^-p*sm-j-l p*sh ^\ 

where the starred quantities are the corresponding tensors as 
measured on stellar-sized objects. Note that the weight function 
with which P* sh and p* sm are measured should be the same as 
that used for the respective object. 

The basic assumption of weak lensing studies is that the av- 
erage source ellipticities vanish, i.e. (^ (s) ) = 0. We also assume 
that ((P g ) _1 A /S ) = 0, so by averaging eq. Q we obtain: 



g = <(p g ry raso > 



(5) 



Thus the expectation value of the quantity (ps) _1 ^ anlso is the 
reduced shear g at the respective point. 

P g is an almost diagonal tensor with similar elements on 
the diagonal. In fact, in the absence of a weight function and a 
PSF its elements would be: P s u = P g 2 = 2 , P s u = P\ Y = 0. We 
can therefore approximate the tensor P g by a scalar quantity: 



Pf = ^tr(P g ) 



(6) 



It has been shown that while the full P g tensor can overestimate 
the shear, this approximation is more conservati ve and will only 
underestimate the true shear iErben et al.l200ll) . 

We therefore use a variant of eq. l|5} as our estimate of the 
local shear at each galaxy's position: 



(7) 



P s 

1 .V 



At this point, we reject those objects from the catalog that have 
a radius equal to or smaller than the stellar locus, are saturated 
stars, or have a final ellipticity of |s| > 1 . 

2.3.8. Rotation of ellipticities 

The ellipticities were measured with respect to the x-axis of 
each image, which for all images runs approximately along the 
east axis. However, the lensing analyses are done in sky coor- 
dinates. We therefore need to transform the ellipticity measure- 
ments so that their position angle is measured relative to the 
right ascension axis. This is done with the transformation 

where <p is the angle of rotation of the image. 

After this step, the 18 single-image catalogs of the HST 
data are merged into one catalog. 

2.3.9. Weighting 

To describe the reliability of its shape measurement, we want 
to assign a weight to each galaxy, based on its noise proper- 
ties. Since our shear estimates are gained from averages over 
ellipticities, a good weight estimate is 



Wi oc 



1 



(8) 
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where cr, is the (two-dimensional) dispersion of image elliptic- 
ities in the ensemble over which is averaged. For intrinsic el- 
lipticities, cr B ~ 0.3 - 0.4 (measured from data). Measurement 
errors cau se cr, to be h i gher t han this. We follow the weighting 
scheme ofEr ben et alllboOll) . and select as the ensemble to av- 
erage over the 20 closest neighbors of the respective galaxy in 
the parameter space spanned by half-light radius r g and magni- 
tude m. 

2.3.10. Final cuts 

After the weighting, we remove objects with |e| > 0.8. Such 
objects would dominate the shear signal, but these are also 
the objects that are most afflicted by noise in the P g ten- 
sor. Additionally, we use only objects for which P g > 0.3. 
This leaves about 2000 objects in both catalogs, which cor- 
responds to 20 galaxies/arcmin 2 for the /-band image and 65 
galaxies/arcmin 2 for the HST image. 



where (f> is the polar angle of the vector The weight function 
Q is determined in terms of w. 

Eq. © is intuitively clear: a lens most often deforms im- 
ages so they align tangentially to the center of mass. An aver- 
age over the tangential components of galaxy ellipticities must 
therefore be a measure of the surface mass. With this inter- 
pretation, M ap is a useful quantity in its own right even if the 
weak lensing assumption, y = (e), breaks down or if part of 
the weight function lies outside the field. 

The imaginary shear component is the cross component: 

T x(0o; 0) = -Im (y(0 Q + #)e- 21 ^) . (11) 

Substituting it for the tangential shear in eq. (|9j yields M* p 
whose expectation value vanishes. Evaluating it analogous to 
M ap can thus be used as a method to check the quality of the 
dataset. 



3. Weak lensing methods 

Weak lensing analyses are based on using estimates of the local 
shear y to reconstruct information on the convergence k, which 
is a dimensionless measure of the surface mass density. In the 
weak lensing regime, k <k 1, so that (e) = g = y/(l - k) « y. 

3.1. Mass reconstruction 

Both the shear y and the convergence k are linear combina- 
tions of second derivatives of the lensing potential, so that it 
is possible to expr ess k as an integral oy er y via the Kaiser- 
Squires Inversion ( Kaiser & Sauires 1993). This method is usu- 
ally not applied directly, as the shot noise introduced by sum- 
ming over individual galaxies (shear measurements) produces 
infinite noise. This can be avoided by first smoothing the shear 
measurements; however, such a smoothing scale introduces 
correlated errors. Another proble m arises from the limited 
field-of-view of any observations. ISeitz & Schneider! J200ll) 
express this as a von Neumann boundary problem, leading to 
the so-called finite-field inversion. We rely on this method for 
mass reconstructions throughout the paper. 

3.2. Mass-aperture statistics 

The m ass-aperture, or M ap , Statistics, developed bv lSchneideJ 
Jl996l) . provides a method with defined noise properties to 
identify mass concentrations. It is based upon the relation 

M ap (0 o )=r d 2 M0o + 0)w(|0|)= f d 2 # rt (0 o ;#)e(|#|).(9) 
Jr 2 Jr 2 

The aperture mass M ap (0o) presents a measure of the average 
convergence, mulitplied with a filter function w, around a po- 
sition 6q in the lens plane. If w(\0\) is a compensated filter, M ap 
avoids the mass sheet degeneracy. The right side of eq. ex- 
presses M ap in terms of the tangential shear y t at the position 
0o + # with respect to 0o: 

yt(0 O ;&) = -<Re (y(0 o + d)^) , (10) 



3.2.1 . Application to real data 

- In order to apply the weight function Q to finite data fields, 
a cut-off radius 6 out should be used, beyond which the fil- 
ter function vanishes. Otherwise, the area of integration is 
not well sampled by galaxy images. A compensated filter 
wQffli), for which w(\&\) = for \&\ > 9 ou i yields a weight 
function Q(\it\) which vanishes beyond the same cut-off ra- 
dius. We use fi l ter an d weight function as introduced in 
ISchneider etail i ll 998b . with I = 1. 



1 




-0.4 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

x = 9 I e out 



Fig. 3. The filter function w (solid line) we used and its cor- 
responding weight function Q (dashed line) shown in units of 

$out- 



- We can sample the shear field only at discrete points, 
namely those where there is a measured galaxy image. In 
the weak lensing regime, the image ellipticity is on average 
a direct measure of the shear, so that we can use the tan- 
gential ellipticity e t as an estimate for the tangential shear 
y t (e t defined analogously to y t ). 
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We thus estimate M ap by using 
nO 2 

M ap (0 o ; flout) = J] em) Q(\0, - 8 \) 



(12) 



where ; runs over all N galaxies within a radius fl ou t from the 
point Oq . 

With the weighting scheme introduced in eq. (|8j, this be- 
comes: 

i, ,a a ^ * Z/gt(g/)trr 2 6(lft-«Dl) 

Af ap (0o; flout) = 71%, . (13) 

Li a~. 



3.2.2. Significances 

Any M ap value is incomplete without an estimate of its signifi- 
cance, i.e. how it compares to the typical noise level of the M ap 
estimator. The signal-to-noise ratio is given by 



S_ _ Map 

N ~ o- M „ 



M m 



^j{Ml v ) - <M ap ) 2 



(14) 



where the denominator should represent the case of no lens- 
ing. The expectation value (M ap ) of M ap then vanishes since 
the galaxies are oriented completely randomly and thus (s t ), 
the expectation value of tangential ellipticities, vanishes. In the 
case of a non-weighted M ap estimator, we obtain: 



(15) 



In the case of no lensing, (e? ) is the one-dimensional variance 
of the intrinsic ellipticities: 



<4> = ^ 

and the signal-to-noise ratio becomes: 

s V2 1,^)2(1^-001) 



^ ^Zi QHWi - O \) 



(16) 



(17) 



For the weighted estimator, one faces the problem that the 
weights 1 /of are in general not completely independent of the 
tangential ellipticity s,j. We assign weights to galaxy images 
by considering the variance of ellipticities of an ensemble of 
galaxy images with similar noise properties. But in general, 
large ellipticities are often noise-afflicted, so that they are as- 
signed a lower weight. The expression for (M^ p ) then cannot 
be simplified ad hoc: 



, , I Zij stjstj o- , 2 cr : 2 QiQ, 



^ijcrfo-f 



(18) 



The significance of a detection is related to the probabil- 
ity that the observed alignment of tangential ellipticities can 
be mimicked by a random distribution of galaxy ellipticities. 
A commonly used possibility to determine the significance is 
therefore to randomize the position angles of the galaxy images 
and calculate M ap of these. This is repeated N mn d times. 



Table 1. Summary of the differences of t he obje ct catalogs used 
for the lensing analyses of Erbe n et alJ (EoOO) and this work. 
The brightness key refers to the magnitude bins used to split 
the sample (note that the magnitudes are not calibrated). 





Erben et al. ("20001 


this work 


masking 




V 


source extraction 


6px > l.Ocr 


3px > l.Ocr 


weighting 




V 


M 


< 1.0 


<0.8 


brightness 


m > 23 


all 






m > 23.67 






22.54 <m< 23.67 






m < 22.54 



With this in mind, we reconsider eq. Jl 8b . Such randomiza- 
tions represent a possible ensemble average, which we denote 
by (. . For each realization, the ellipticity modulus remains 
the same, only the orientation changes. In this case, the weights 
also remain the same and we can simplify the expression: 



since, as the ellipticity modulus is fixed, 

/ \ _ N 2 * 

\£u'£tj/0 - — ^— Oij . 

The signal-to-noise ratio for the weighted estimator is there- 
fore: 



s_ = ./o y em ^ Qil6i ~ ° ol) 
N ' ^i\si\ 2 o-T 2 Q 2 (\0i-9o\) 



(19) 



To check the validity of the assumptions we made, we com- 
pared results from randomizations and this analytic formula 
and found them to be equivalent. 

4. Re-analysis of the CFHT data 

In TableU] we summarize the differences between the analyses 
of lErbenetaH l2000) and this work. 

4.1. Mass reconstruction 

In order to apply the finite-field inversion (Sect. l3.ll . we further 
cut the /-band image to avoid the ripple-like reflection artefact 
at the eastern edge altogether. This narrows the available field, 
but avoids problems at the boundaries. A resulting mass recon- 
struction is shown in Fig.0] 

Abell 1942 shows up prominently in the top half of the 
field, with the peak of the mass map centered approximately on 
the cD galaxy. In the lower half of the image, there is a second, 
albeit lower peak at the same position as detected in tErben et alJ 
( 2000). Relative to the peak k of A1942, our dark clump signal 
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Fig. 4. M ass reconstruction of the CFHT /-band image ac- 
cording to lSeitz & SchneideJll200l|) . Solid (dashed) lines give 
positive (negative) k contours, starting at ±0.02 and in-(de- 
)creasing in 0.01 intervals. The smoothing scale is 60". 



is slightly larger than given bv lErben et aD J200(A . However, 
north of the dark clump, there is a "hole" - a region of signif- 
icantly negative k values. Although this is a somewhat discon- 
certing result, it must be stressed that k is underestimated in the 
whole field due to the mass sheet degeneracy and the cluster in 
the field (two of the three field boundaries close to the cluster 
and well within its extent display nearly vanishin g k va lues). 
Unfortunately, the original analysis of Erb en et al.l 112000) only 
investigates regions of positive k, so that this result cannot be 
compared. 

4.2. M ap analysis 
4.2.1. Complete sample 

We perform the M ap analysis as described in Sect. l3.2l and eval- 
uate the M ap statistic on a grid with grid spacings of 2". The 
result is shown in Fig. EI The dark clump signal is seen signifi- 
cantly at all filter scales, but it is particularly strong for the 120" 
filter, where it reaches a peak significance of 5cr. In the o ther 
filters, the significance is at the 3.5<x level, as in lErben et alJ 
(200B- 



Abell 1942 is detected only wea kly for large fi lter scales, 
a result that is consistent with Erben et al.l (I2000I) . The mass 
reconstruction from the previous paragraph demonstrates why 
the dark clump reaches a much higher M ap significance than 
the cluster: at the dark clump position, the negative part of the 
filter function (Fig. [3} is evaluated largely at the position of the 
hole, thereby boosting the signal. The same in reverse is true 
for the hole itself: its significance is boosted by its proximity to 
the Dark Clump. Yet its significance remains lower than that of 
the dark clump. 

4.2.2. Rough redshift dependence 

On average, the more distant a galaxy, the fainter it is. By in- 
troducing magnitude cuts we split the galaxy sample into three 
parts of about 660 galaxies each with different mean redshift. 
This is a very crude redshift distinction, but should reveal any 
trend of lens strength with redshift. 

The results of this analysis are shown in Fig. |6] We see 
that the dark clump signal stems mostly from faint galaxies, 
which supports the notion that this is a high-redshift object. 
However, these are also those objects that are most subject to 
noise effects. 

At the 120" filter scale, there is also a 3<x contribution 
from bright galaxies. Assuming that the lensing mass is indeed 
a high redshift object, these bright galaxies are unlikely to be 
at higher redshifts. Thus, this is probably not a lensing signal. 
Yet, this "contamination" can explain the high signal-to-noise 
ratio we see at this filter scale. 

For Abell 1942, there is a strong signal at the smallest fil- 
ter scale, centered on the cD galaxy which exhibits a strong 
lensing arc. It may well be that at these radii we are not in the 
weak lensing regime any more and the tangential alignment is 
already rather distinct. In the other filter radii, there is no par- 
ticularly strong signal. This might be due to the generic weight 
function which is not adapted to the NFW profile. 

Our plot does not show the negative M ap contours to avoid 
overcrowding. Unlike to the dark clump, all three magnitude 
bins contribute to the "hole". 

4.2.3. M ap cross component 

In Sect. 13.21 we argued that M ap , i.e. M ap calculated with the 
cross component of the shear instead of the tangential one, 
must vanish. By checking the validity of this assumption in the 
dataset, we can identify possible problems. 

The results are shown in Fig.0 Particularly at small filter 
scales, there are some positive and negative peaks with signifi- 
cances >|3cr| (negative values are indicated by negative signal- 
to-noise ratios). However, most of the peaks are at the edge 
of the field, where a part of the weight function lies outside 
the field. While (the real part of) M ap retains its justification at 
these places as being simply a measure of the tangential align- 
ment, we can no longer assume that M ap vanishes. 

There is no peak with a significance larger than 2<x in the 
vicinity of the dark clump, so the detection passes this test well. 
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Fig. 5. Results of a weighted M ap analysis of the /-band data. Shown are positive (negative) S/N contours as solid (dashed) lines 
starting at + 2. Oct in 0.5<r intervals. The filter scales from left to right are 80", 120", 160", and 200". The dark clump signal 
shows up strongest in the 120" filter, with a peak significance of 4.9<x. 
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Fig. 6. The same as the previous figure (Fig. [5}, but with the galaxies divided into three magnitude bins and only positive 
contours. White contours correspond to the brightest galaxies, gray contours to those of medium brightness, and black contours 
to the faintest. 



4.3. Radial profile pendence of the m ean tangential ellipticity shown in Fig. 9 of 

lErben et alJl200(]ll . 

Despite se veral differences to the lensing analysis of We determine the position of the dark clump from the M ap 

Erb en et alJ (EoOO). the M ap results agree at least qualitatively, peak at the 120" filter scale, where the signal is strongest. We 
We are curious whether we can also reproduce the radial de- find it to be a = 14 h 38 m 21.6 s , 5 = 3° 32' 43'.' , which is at a 
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Fig. 7. Results of a M* p analysis of the /-band image. Shown are both positive (white) and negative (black) contours, starting at 
±2<x in 0.5cr increments, for filter scales of 80", 100", 160", and 200" (left to right). 
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Fig. 8. The mean tangential ellipticity (e t ) relative to the posi- 
tion of the Dark Clump as a function of distance from it (from 
the ground-based data). Each point corresponds to the mean 
of a 20" wide bin, where filled circles denote the complete 
sample, crosses the faint galaxies, triangles the medium bright 
ones, and squares the bright galaxies. The error bars represent 
the standard deviation of (e t ) estimated from randomization of 
the position angles of the galaxies. Note that the first bin is very 
sensitive to the choice of centroid - its deviation from the shear 
profile in the latter bins is thus not problematic. 



distance of 18'.'6 from the position given bv Er ben et aPfcOOOl) . 
and thus just at the lcr level they give for the uncertainty of the 
centroid's position. 

The tangential ellipticity relative to this position is calcu- 
lated for each galaxy within 200". They are then binned ac- 
cording to their distance from the dark clump and the weighted 



mean is calculated for each bin. To estimate the standard de- 
viation, we randomize the position angles of these galaxies 
1000 times and calculate the mean tangential ellipticity each 
time, thus gaining an estimate for the standard deviation. This 
analysis is done for the complete galaxy sample as well as the 
three samples split according to brightness. The results of it are 
shown in Fig. [8] 

Particularly for the faintest galaxies, we find positive values 
out to 120". This agrees well with the strong shear signal 
seen for these galaxies. We can also identify the cause of the 
signal seen for bright galaxies at the 120" filter scale: the two 
significantly positive bins at 90" and 110" (the filter function 
employed assigns the highest weight around a radius 8 out / V2). 
For the medium bright galaxies (e t ) is largely consistent with 
zero. 

Compared to Erbe n et al who measure (e t ) ~ 0.06 

at 100", we find a higher value (» 0.1). On the other hand, 
we find positive values only out to 120" rather than 150". And 
since the centroid positions do not coincide, the inner two bins 
are not comparable. Yet, we can also be confident that the sig- 
nal is not caused just by a few galaxies. 

4.4. Summary 

We have successfully confirmed the weak lensing signal seen 
in two sets of CFHT observations (our re-analys is of the V- 
band data are not shown here but agree well with Erben et all 
(2000)). We show that the alignment signal comes from faint 
galaxies, which supports the hypothesis that it is caused by a 
lensing mass at high redshifts. One must keep in mind, how- 
ever, that these are also those objects most affected by noise. 

With several variations of the catalog that enters the M ap 
analysis, we tested that the detection of the dark clump is resis- 
tant against these and consistently recovered at all filter scales. 
It reaches a peak significance of about 5cr, although this sig- 
nal is contaminated by a tangential alignment of bright objects, 
which is unlikely to be a lensing effect. 
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Table 2. Overview of the M ap values and signal-to-noise ratios of the peaks found in the analyses of the CFHT data. We also 
give the number of galaxies located in the aperture (AO, and the (weighted) average galaxy ellipticity dispersion cr e , as well as 
the offset to the assumed centroid position, measured in the sky coordinate system. In two cases, there was no M ap peak in the 
vicinity of the dark clump - we then quote the values at the reference position (in italics). 







N 




Aa["] 


M["] 




SNR 


|e| < 0.8, 


80" 


113 


0.34 


28 


16 


0.096 


3.7 


weighted, 


120" 


268 


0.34 








0.078 


5.1 


all galaxies 


160" 


468 


0.35 


-16 


-2 


0.046 


3.7 




200" 


648 


0.34 


58 


-12 


0.034 


3.5 


|e| < 0.8, 


80" 


54 


0.41 


22 


26 


0.120 


2.7 


weighted, 


120" 


110 


0.39 


4 


4 


0.109 


4.1 


m > 23.67 


160" 


167 


0.38 


36 


18 


0.092 


4.3 




200" 


241 


0.39 


40 


-22 


0.088 


4.5 


|e| < 0.8, 


80" 


33 


0.34 


40 


8 


0.118 


2.6 


weighted, 


120" 


86 


0.35 


6 


-18 


0.058 


2.0 


22.54 < m < 23.67 


160" 


159 


0.36 


-26 


2 


0.031 


1.5 




200" 


222 


0.34 






-0.003 


-0.1 


|e| < 0.8, 


80" 


35 


0.29 


38 


4 


0.084 


2.2 


weighted, 


120" 


74 


0.29 


16 





0.075 


3.5 


m < 22.54 


160" 


126 


0.30 


-14 


-20 


0.046 


2.3 




200" 


170 


0.30 


-48 


34 


0.026 


1.6 



5. Analysis of the HST data 

The HST catalog extends to fainter and thus more distant galax- 
ies than the ground-based catalog. If the alignment signal found 
in the ground-based data does indeed stem from a mass con- 
centration, its lens strength should increase with source red- 
shift. Additionally, since the HST is space-based and thus not 
afflicted by atmospheric seeing, its ellipticity measurements are 
more reliable than those from ground-based telescopes. If the 
ground-based signal is not just a noise peak, the signal should 
therefore be even stronger in the HST images. 

5.1. Mass reconstruction 

Fig. |9] shows the results of a mass reconstruction of the in- 
ner rectangle of the HST data field. The dark clump shows 
up prominently, slightly westward of the position found in the 
ground-based analyisis. 

5.2. M ap analysis 

The results of the weighted M ap analysis of the complete cat- 
alog is shown in Fig. ^| and summarized in Table [3] Indeed, 
we find a peak at approximately the same position as in the 
ground-based images in the 120" filter, but with only 2.9<x sig- 
nificance it is considerably weaker. At the smallest filter scale 
used (80"), this peak diminishes and dissolves into two peaks. 

Just as for the ground-based data, we split the catalog into 
three bins according to brightness to probe the redshift de- 
pendence of the lens strength. With cuts at m = 25.46 and 
m = 26.43, the three samples contain approximately equal 




Fig. 9. The mass density map of the inner region of the 
HS T image recon s truct ed according to the method given in 
ISeitz & Schneider! J200ll) . Shown are k contours starting at 
+0.02 and in-/decreasing in 0.01 decrements. The smoothing 
length is 30". 



numbers of galaxies. The respective M ap analyses are shown 
in Fig. \^\ The alignment signal is carried by the galaxies in 
the bright and in the faint bin, but there is a lack of signal in 
the medium-bright bin. The galaxies that were detected in the 
ground-based image should be mostly contained in the HST's 
bright bin. Although the M ap values differ by a factor of about 
2, the measurement in the bright bin therefore confirms that 
there is tangential alignment around the dark clump candidate. 




Fig. 11. The same as the previous figure (Fig.llOI. but with the galaxies divided into three bins. White contours correspond to the 
brightest galaxies, gray contours to those of medium brightness, and black contours to the faintest. The contours start at 1 .5cr and 
increase in 0.5cr increments. 
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However, the lack of alignment in the medium bin is difficult 
to explain with the lensing hypothesis. 

5.3. M ap cross component 



0.4 



As described in Appendix lA.il it is difficult to judge the qual- 
ity of the anisotropy correction of the HST data. A faulty 
anisotropy correction is likely to cause a non-vanishing M* p 
component. Also the CTE problem of the camera might do so. 
We therefore calculate M ap with the shear cross component in- 
stead of the tangential component; the result of this analysis is 
shown in Fig. 1121 

Indeed, there is a 3<x peak roughly coincident with the 
group of galaxies close to the dark clump, and one close to 
the border of the image. For the latter one, a large part of the 
aperture is evaluated outside of the field, so that M ap is not nec- 
essarily expected to vanish in this case. This argument applies 
only weakly to the first one, which is affected only by masks of 
stars. 

The stellar anisotropy increases with distance from the cen- 
ter of the chips, cf. Fig. IA.4l At the edges of the chips we also 
expect the largest deviations of the fit from the actual values. 
Therefore, the quality of the anisotropy correction should de- 
crease with increasing distance from the chip center. To test this 
effect, we reject those objects which are more than 700 pixels 
from the center of their chip (343 objects) and redo the M ap 
analysis. The result differs only very little from the previous 
one. 

We also calculate the M ap cross component for the three 
brightness bins, Fig. El We find that the M* p peaks are caused 
solely by the bright galaxies. If they were caused by an insuf- 
ficient anisotropy correction, we would expect these peaks to 
show up for all three bins. Additionally, as the brightest ob- 
jects are generally also the largest, the effect of the anisotropy 
correction is smallest for these. 

These results indicate that the M ap peaks are unlikely to be 
caused by a poor anisotropy correction. Rather, they seem to be 
caused by an intrinsic alignment of some of the bright galaxies. 

5.4. Radial profile 

Just as described in Sect l4.3l we calculate the average tangen- 
tial ellipticity (s t ) as a function of the distance from the Dark 
Clump position. The position of the strongest lensing signal is 
a = 14 h 38 m 20.3 s , 5 = 3° 33' 10", as measured from the peak 
in the 120" filter of the complete sample. This is 33'.'3 from the 
position we measured in the ground-b ased data an d 5 l'.'O from 
the position originally cited by Er ben et alJ (12000). 

The radial profile is shown in Fig. El F° r the complete 
sample of galaxies, (e t ) is positive or very close to zero between 
20" and 140" radii. For this range, it is inconsistent with zero 
only for radii between 60" and 100". 

For bright galaxies, there is some excess (e t ) between 20" 
and 100", which is clearly the cause of the peak we find. For 
the medium bright galaxies, (e t ) is largely consistent with zero 
within in the lcr error bars. However, the value is positive only 
in two bins. For the faintest galaxies, <e t ) is positive over a 




Fig. 14. The mean tangential ellipticity (s t ) to the position of 
the dark clump (determined from the M ap analysis of the HST 
data) as a function of distance from it. Each point corresponds 
to the mean of a 20" bin, where filled circles denote the com- 
plete sample, crosses the faint galaxies, triangles the medium 
bright ones, and squares the bright galaxies. The error bars rep- 
resent the standard deviation of (e t ) estimated from randomiza- 
tion of the position angles of the galaxies. 



fairly large range of radii, but with varying significance and no 
clear resemblance of a shear profile. 

We also calculate the radial profile around the dark clump 
center we found in the ground-based images; it is shown in 
Fig-El F° r the complete galaxy sample, it is largely consistent 
with zero but with a trend to positive values. But for the bright 
galaxies, (e t ) follows a typical shear profi le between 20" and 
100", with (e t > * 0.025 at r * 100" lErben etai] d2000h use 
the radial profile, namely (e t ) m 0.06 at r — 100", to deduce 
a mass estimate of the Dark Clump. This further illustrates the 
difference of a factor of 2 between the signal strengths in the 
ground-based and space-based image, which relates directly to 
the mass estimate. 

However, for the medium bright and faint galaxies, there is 
no obvious NFW-like trend in the profile. 

5.5. Projected galaxy density 

Neither lErben et alJ fcOOOt /-band) nor ldravet all fcOOll H- 
band) find a significant overdensity of galaxies which could be 
associated with the Dark Clump. With the substantially deeper 
HST data, we can probe the galaxy number density to fainter 
magnitudes. In Fig. El we present surface number densities 
for the complete galaxy sample, as well as split up according 
to brightness, as done before for the lensing analysis. 
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Fig. 12. The M* p analysis of the HST data. Shown are both positive (solid) and negative (dashed) contours, starting at + 1 .5cr in 
0.5<x increments, for filter scales of 80", 100", 120", and 140" (upper left to lower right). 
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Fig. 15. Same as Fig. [H] but with the distance measured rel- 
ative to the position of the dark clump (as measured in the 
ground-based data). 



The bin of brightest HST galaxies corresponds roughly to 
those objects detected also in the ground-based image. It is 
quite puzzling that rather than an overdensity of galaxies at 



the dark clump, the galaxies seem to form a ring around that 
position. However, this is hardly significant. 

The most significant feature we find is a 5<x galaxy over- 
density in the bin of medium bright galaxies, located about 1' 
from the Dark Clump. Given that these objects are fainter than 
the ones that carry the lensing signal in the ground-based data, 
this is rather unlikely to be associated with a lensing mass. 
However, it could explain the alignment signal seen in the 
faintest HST bin. It would be highly interesting to investigate 
whether this overdensity is present also in color space; unfor- 
tunately, the currently available data sets are not deep enough. 

5.6. Summary 

We have shown that we also detect tangential alignment around 
the Dark Clump in the HST data. However, it is considerably 
less significant than in the ground-based data. It is particularly 
intriguing that we do not detect alignment in the medium bright 
HST galaxies. 

We refrain from trying to deduce a mass for the dark clump 
from these measurements, as it is obvious that our results are 
not unambiguous. As an upper limit, the HST data suggest that 
the mass estimate of the dark clump has to be corrected at least 
by a factor of 1/2 compared to the original value. 

A comparison on an object-to-object basis may be able to 
shed some light on the question whether the ground-based lens- 
ing signal is only a noise peak or whether there are other, un- 
corrected problems with the HST mosaic. Such a comparison 
is carried out in App.lBl 




Fig. 13. The same as Fig.^] but with the galaxies split into three samples according to brightness. The color coding is identical 
to that of Fig. ^2 The contours start at ±1.5cr and in-/decrease in 0.5<x increments. 

Table 3. Overview of the M ap values and signal-to-noise ratios of the peaks found in the analyses of the HST data. We also give 
the number of galaxies located in the aperture (AO, and the (weighted) average galaxy ellipticity <x £ , as well as the offset to the 
assumed centroid position, measured in the sky coordinate system. For the small filter scales (60" and 80"), there were often 
more than one peak in the vicinity of the dark clump position. For the medium bright bin, there is no peak close to the dark clump 
in the largest two filter scales - we then quote the values at the reference position (in italics). 
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Fig. 16. Projected number density of galaxies in the HST image. The upper left image shows contours of surface densities, 
smoothed with a Gaussian of width <x = 15", for all galaxies in the catalog, starting at 70 arcmkr 2 , and increasing in 10 arcmkr 2 
intervals. For the other three plots, the galaxy catalog was split as before into magnitude bins. Bright galaxies are shown in the 
upper right panel, medium bright ones in the lower left, faint ones in the lower right. The contours for these panels start at 
25 arcmin~ 2 , and increase in 5 arcmkr 2 steps. 



6. X-ray analysis 

If the dark clump is a virialized massive dark matter concentra- 
tion, one would expect that the baryons trapped by the poten- 
tial well and not in galaxies would have been shock-heated in 
the gravitational c ollapse and therefore would emit in X-rays. 
Erb en et alJ (EoOO) present the analysis of the ROSAT HRI im- 
age of the A 1942 field. They found a hardly significant X-ray 
detection at a position 60" away from the weak lensing mass 
peak position reported in that paper. 

Since then, the field of A 1942 was imaged by the Chandra 
X-ray observatory in March 2002 (PI: Garmire). It was ob- 
served for 58 ks with the ACIS-I detector configuration. We 
retrieved the observation from the archive and analysed it using 
standard techniques with the CIAO v3.1 package. 1 and recipes 
from the "AOS Recipes: Clean the Data" web page. 2 We use 
the standard set of event grades (0,2,3,4,6) and restrict our anal- 
ysis to the 0.5-8.0 keV energy band. We detect sources with 
CIAO's WAVDETECT routine using wavelet scales 1, 2, 4, 8 
and 16. Only two sources are detected within 100" of the CFHT 
data lensing M ap centroid. These two sources are found to be 
point sources with clear optical counterparts in the optical data. 



1 Chandra Interactive Analysis of Observations (CIAO), 
http://cxc.harvard.edu/ciao/ 

2 http://www.astro.psu.edu/xray/acis/recipes/clean.html 



Visual inspection does not reveal any obvious extended 
source at the position of the dark clump. Nevertheless, in or- 
der to check if there really is any extended source, we cre- 
ated a diffuse emission image from our original image excis- 
ing the data where WAVDETECT detected sources and fill- 
ing the holes according to CIAO's threads. We also created 
a "blank-sky" background image to better estimate the back- 
ground contribution at the position of the dark clump. We 
smoothed the point-source-removed image with and without 
background subtraction with various algorithms and scales but 
no obvious extended source was detected. The most significant 
source consistent with the broad position of the dark clump is 
a 1.5cr source (in a circular aperture of radius 45") whose co- 
ordinates are a = 14 h 38 m 22.7 s , 6 = 3° 32' 57". We measure 
a count-rate of (5.9 ± 4.0) x 10 4 counts/s in the 0.5-8.0 keV 
energy range in a circular aperture of radius 45" which corre- 
sponds to a unabsorbed flux of (5.3 ± 3.4) x 10" 15 erg crrT 2 s _1 
assuming an incident spectrum of T = 3keV and a local hydro- 
gen column density of Nh - 2.61 x 10 20 crrT 2 . We also measure 
the count rate in various apertures centered at the positions of 
the CFHT and HST data lensing centroids. The measurements 
are of lower significance compa red to the previous one. We 
have fitted a standard beta profile (Cavalie re & Fusco-Femianol 
119781) to the azimuthally average radial profile centered at the 
position of this potential source. We obtain best values for the 
core radius and beta parameter (slope decline at large radius) 
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Fig. 17. A part of the Chandra image overlaid on the HST image. The point sources were removed from the image, and 
it was smoothed with a Gaussian of 5". The contour levels start at 3.6 x 1CT 14 erg cirT 2 s _1 and increase in steps of 
0.8 X 1CT 14 erg cm -2 s _1 . An increase in X-ray emission due to Abell 1942 (north of the dark clump) is visible. In the vicinity of 
the dark clump, there is a 1.5cr source at a = 14 h 38 m 22.7 s , 5 = 3° 32' 57" (marked with x) in a circular aperture of radius 45" 
(indicated). 



of 5" and 0.55, respectively. These values indicate that if there 
really is an extended source, it is rather compact with a small 
core radius. Although, one has to keep in mind that these best 
fit parameters are highly uncertain due to the low number of 
counts. The total count rate predicted integrating the best-fit 
mo del would be 9.4 x 10~ 4 counts/s. 

lErben et ail 12000) reported a hardly significant (3.2<x) de- 
tection of X-ray emission in the dark clump area in the ROS AT 
HRI image of this field. We have measured the flux in the 
Chandra image at the same position, which is very close (15" 
away) to the previously reported source. We find an unabsorbed 
flux of (6.8 ± 3.6) x 10~ 15 erg cm~ 2 s _1 in a 45" aperture. This 
value is slightly higher than the previous one only due to the 
inclusion in the aperture of another faint source situated at the 
NE and not included in the mesurement at the previous posi- 
tion. Therefore, we do not confirm the flux measurement of the 
ROSAT HRI source. 

Given the faintness and measured uncertainties of our pos- 
sible detection, there is little point in speculating about the lu- 
minosity and gas mass of our possible detection. Even if this 



source was real and its real flux was the highest allowed by the 
data, it would not have enough mass to be considered a rich 
cluster by its X-ray properties. Overall, the Chandra image of 
this source indicates that if there is a dark matter concentration 
producing the lensing signal, this concentration does not con- 
tain the expected mass of hot gas that would be expected for its 
lensing signal if it were similar to the clusters of galaxies we 
have observed so far. 

7. Conclusion 

Tangential alignment around the dark clump was detected in 
three datasets, which differ in the filter, camera, and telescope 
used. It can therefore be ruled out that the alignment is caused 
by instrument-specific systematics. However, the significance 
of the detection is lowest in the space-based dataset. If the 
alignment were due to lensing by a matter concentration, the 
highest alignment signal would be expected in the HST data, 
as it extends to more distant galaxies which should be more 
affected by the distortions due to lensing. 
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The significant alignment signal in the ground-based data 
is caused mainly by faint galaxies, for which the individual 
shape measurement are uncertain due to background noise. We 
show in Add. 151 that the shape measurements agree very well 
on average when comparing ground-based and space-based 
data. Thus, also the amplitude of an alignment signal should 
be comparable. Considering the high alignment signal in the 
two ground-based datasets, it seems unlikely that background 
noise boosted the signal in both cases. Yet, for the /-band im- 
age, there is some indication of such a "conspiracy", as the av- 
erage tangential ellipticity toward the dark clump of the objects 
compared is significantly larger for the ground-based measure- 
ments (see lB.4> . 

There are several issues involved which actually make 
weak lensing analyses more difficult to apply to HST data 
(small field-of-view, complicated PSF structure). But various 
tests (App.lAland lB.4> have not revealed any bias these prob- 
lems might cause to the ellipticity measurements. The weaker 
alignment signal in the HST data is therefore not a result of 
these systematics. 

For the original detection claim, further evidence for a 
lensing mass was lent by the detection of X-ray emission in 
a ROSAT HRI image. However, follow-up observations with 
Chandra measure only a tenth of the flux originally measured 
in the ROSAT data, making this more likely to be a spurious 
detection. Thus, if it is a lensing mass, the dark clump would 
be truly "dark", not just in the optical. 

We cannot give a definite answer to th e question of the na - 
ture of the Dark Clump candidate found bv lErben et alJd2000l) . 
It has neither been possible to prove that it is a Dark Matter 
halo, nor that it is a statistical fluke, caused by noise in either 
the intrinsic galaxy ellipticities or the measurement process. 
Our analyses assert that there is indeed tangential alignment 
around the dark clump. However, it remains unclear whether 
this is caused by a lensing mass, is just a chance alignment, or 
a combination of the two. 

The lack of signal in the medium magnitude bin of the HST 
data and the lack of X-ray emission yield the first hypo thesis 
more unlikely (at least compared to lErben et al. (2000)), but 
not impossible if one is willing to accept the existence of very 
dark, massive halos. 
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Appendix A: Testing for systematics in the HST 
data 

Despite the lack of seeing due to the Earth's atmosphere, weak 
lensing analyses of HST data are not applied straightforwardly. 
Due to the small field-of-view, there are only few stars available 



for the anisotropy correction; the pixels of the WFPC2 cam- 
era undersample the PSF, and the camera has a notable charge 
transfer inefficiency. In the following, we investigate the effects 
and possible bias from these problems. 

A. 1. Anisotropy correction 

iHoekstra et alJ dl998l) showed that a better anisotropy correc- 
tion (eq. 0) of WFPC2 images can be achieved if the weight 
function applied to the stellar images is adapted to the size of 
the object to be corrected. 

However, in practice, the ellipticity measurements become 
increasingly noisy for larger radii of the weight function. Along 
with the low number of stars, this causes fits to the anisotropy 
kernel to be both noisy and ill-described by a simple polyno- 
mial. 

We tested three possible methods of correcting for the 
anisotropy as measured from the stars: 

1 . the anisotropy kernel q* fit to measurements with a weight 
function of stellar size; 

2. the anisotropy kernel q*(r) fit to measurements with a 
weight function of the objects' size r; 

3. a factor a is applied to q* to minimize the difference to the 
stellar ellipticities measured at the objects's size. 

The third method acknowledges that the anisotropy kernel is a 
function of object size, but avoids the uncertainties from fitting 
a polynomial to the noisy measurements with a larger weight 
function. However, it can only compensate for variations in the 
PSF that are uniform across the chip and of the same magnitude 
in both ellipticity components. 

The anisotropy patterns produced by these methods are in- 
deed notably different. But the resulting shear estimates s differ 
only on a percent level. Accordingly, the M ap measurements 
change only minutely with a different anisotropy correction. 
We used the third method as a compromise solution. 

The fits to q* are illustrated in Figs. IA.ll - lA~3l For each 
chip, about 30 stars are available for the anisotropy correction. 
This is a very small number, and as can be seen from the 
Figures, they are not evenly distributed over the chip. It 
is therefore difficult to judge the quality of the anisotropy 
correction. 

Alternatives? It is common to use globular cluster fields 
for the anisotropy correction rather than the few stars contained 
in the field itself. With the large number of stars, the anisotropy 
can be very well determined. However, such fields have to be 
chosen carefully: they should have been taken temporally close 
to our field, have a similar dither pattern, and should have been 
taken in the same filter. As it turned out, no images of star clus- 
ters, M31, M33, the Magellanic Clouds, or galactic fields were 
taken with the F702W filter in 2000 - 2002. 

Instead, we retrieved images of a field within the globular 
cluster 47 Tuc (NGC 104) taken in the F606W filter on Oct. 
19th, 2002. We found the anisotropy patterns of these images 
to be similar to the ones in our images, but with some no- 
table differences. Using these patterns for the anisotropy cor- 
rection did not reduce the ellipticity dispersion of the stars (see 
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Fig. A.l. The anisotropy correction for Chip 2. On the left, the magnitudes of the 850 objects left in the catalog at this stage are 
plotted vs. their half-light radii. 33 objects from the stellar sequence were used to determine the stellar anisotropy kernel. The 
orientation and modulus (multiplied with a factor of 2000 [pixels]) ofx* are shown in the middle top box. The fitted polynomial, 
evaluated at the position of these stars, is shown in the top right box. In the middle lower box, the ellipticity components of these 
stars before the anisotropy correction are shown; in the lower right box they are shown after correction. 



Fig. IA.4> . At least for the centers of the images, the previ- 
ous method is therefore superior. However, due to a lack of 
stars, the anisotropy correction is largely extrapolated toward 
the edges of each chip. These are also the areas with the largest 
anisotropy (as apparent from the figures), so that we are uncer- 
tain of the quality of the anisotropy correction applied to them. 



A. 2. Charge Transfer Efficiency 

WFPC2 has a considerable Charge Transfer Efficiency (CTE) 
problem. Stars bright enough to show diffraction spikes also 
show typical "CTE tails" in the direction opposite to the read- 
out direction. In fainter objects, these tails are not distinguish- 
able by eye. Unfortunately, this effect is difficult to compensate 
for an d has been studied only little for extended objects (Riess 
2000). 

The imperfect CTE of the WFPC2 could potentially bias 
the measured ellipticities, since it results in charge being de- 
pleted in some pixels and added into others. Unlike for the pho- 
tometry of point sources, there is not yet a correction procedure 
to account fo r this effect for shape measurements. According to 
iRiessI J2000l) . for extended objects, the A deficient CTE causes 
charge loss in the pixels closer to the read-out amplifier but 
adds this charge to the pixels farther from the amplifier. To first 
order, this just causes a slight displacement of the galaxy's cen- 
troid. But since depletion and addition cannot be expected to 
be symmetric, it may also distort the galaxy's shape. If so, it 



would bias the si components, as this is measured along rows 
and columns. 

From studies on how CTE affects photometry it is known 
that it increases with distance from the read-out amplifier. This 
is also visible in our images: stars with large y-coordinates have 
more pronounced CTE trails. Also, the charge loss due to CTE 
depends on the background level of the image. A high back- 
ground effectively suppresses CTE losses. The images of our 
dataset have a background corresponding to about 35 electrons 
per pixel, which reduces CTE loss significantly, at least for 
photometry. Lastly, the relative CTE losses are largest for faint 
objects. 

To investigate any possible bias due to CTE, we divide the 
galaxy catalog into four bins according to the original y posi- 
tion, which gives its distance from the amplifier. Unlike for the 
lensing analysis, we are not interested in the ellipticity with re- 
spect to the rectascension axis, but to the original x-axis of the 
image. We denote these with e,-. For each bin, we calculate the 
weighted mean of s\ and si- The results are shown in Fig. IA.5l 

If the CTE would affect galaxies similar to stars, i.e. it 
causes them to trail and thus elongates them in the s\ direc- 
tion, we would expect that {e\) is consistent with zero for the 
first bin and then decreases with increasing distance from the 
amplifier. This is clearly not the case for any of the brightness 
bins. The scatter in {tt\ ) about zero is comparable to that in (£2), 
which should not be affected by CTE. The deviations from zero 
in some bins may well be due to the fact that there is some de- 
gree of tangential alignment in the field, so that (e\) and (§2) 
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Fig. A.3. Same as Fig. IA.ll but for Chip 4. 30 objects out of 770 were used for the fit. 



are not necessarily zero. Yet, Fig. [A3] excludes a notable bias surements from the ground to those from space. Ground-based 
due to CTE. shape determinations rely on an accurate compensation of the 

smearing due to the Earth's atmosphere, while the HST ob- 
servations are hampered by the small image size, and the un- 
dersampling of the PSF In the previous section, we confirmed 
with the HST data the presence of tangential alignment around 
the dark clump candidate, which implies that ellipticity mea- 
surements are to some degree comparable. However, we failed 



Appendix B: Comparison of ground-based and 
space-based measurements 



With the two datasets - the /-band image and the HST im- 
age - we have the opportunity to directly compare shape mea- 
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Fig. A.4. Anisotropy correction of Chip 3. In the top left field, we plot the ellipticities of the 33 stars from the stellar sequence (as 
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to confirm with the HST data the amplitude of the alignment 
signal. 

A direct comparison of the objects common to both cata- 
logs tests the reliability of ground-based shape measurements 
and may also help to identify any systematics in either dataset. 
The ultimate goal of this comparison is to find the cause of the 
discrepancy in the significance of the alignment signal between 
the two datasets. 

B.1. Catalog correlation 

As the HST data was astrometrically calibrated by using a ref- 
erence catalog extracted from the /-band image, objects present 
in both images can be identified by their sky coordinates. 

To correlate the catalogs, we searched for objects within 1" 
of objects detected in the respective other catalog. This radius 
was found to be the optimal balance between a high number of 
matched objects and a low rate of double detections. 

The catalogs used are the same as those for the lensing anal- 
ysis, except that objects with Pf < 0.3 in the ground-based 
catalog are also considered. 

Within the field covered by the HST mosaic, there are 507 
objects in the /-band catalog. We find one HST counterpart for 
350 of these objects, and two or three counterparts for 17 ob- 
jects, as the HST is able to resolve very close objects which 
were identified as single objects in the ground-based data. For 
140 objects, no counterpart was found. For most of these, this 



results from the large areas that were masked in the HST im- 
ages. 

Fig. lB.ll illustrates the matched objects in an r g -m diagram 
of the ground-based data. There is no apparent trend as to for 
which objects we are more likely to find a counterpart in the 
space-based data. In particular, it is not less likely to do so for 
faint objects. This shows that the catalog was only very little 
contaminated by noise detections. Since close objects are often 
unresolved in the ground-based image, objects with two HST 
counterparts have on average a larger radius than those with 
one counterpart. 

B.2. Statistical properties 
B.2.1. Ellipticity measurements 

The basis of weak lensing analyses are the shape measurements 
of faint galaxy images. But the fainter an object is, the more 
difficult the shape determination is. Our dataset provides an 
ideal opportunity to test the reliability of shape measurements 
of ground-based data compared to space-based measurements. 

In Fig. IB.2I we compare the ellipticity measurements in 
the space-based data to those of the ground-based data. To 
thes e points, we fit a linear function (algorithm fitexy from 
IPress et alJll99l to account for errors in both coordinates), 
where we employ the same weighting scheme as applied to the 
lensing analysis. We find: 
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Fig. B.2. Comparison of the ellipticity measurements from the HST images to the /-band image. Objects for which Ae < 0.2 [see 
eq. JB.1H are shown as solid circles, objects with Ae > 0.2 are shown as open circles. The line gives the best linear fit. The error 
bars of individual points are omitted to avoid crowding. 



£i,hst = m\ ■ £1,1 + bi e 2 , hst = m 2 ■ £2,i + b 2 

mi= 0.94 ±0.14 m 2 = 0.97 + 0.15 

bi = 0.017 + 0.027 b 2 = -0.001 ± 0.027 

Although the scatter is rather large and a comparison on an 
object-to-object basis not feasible, the general agreement for 
the is very good, particular for the s 2 component. For the ei 
component, both slope and offset b\ differ from unity by about 
0.5<x. If the HST ellipticity measurements are influenced by the 
CTE problem or by the resampling of the drizzle algorithm, 
we would indeed expect some bias in the e\ component, which 
is measured along the rows and columns of the CCD chip. 

The scatter seen in the correlation of both components is 
of a similar order of magnitude, so that the deviations in the 
ellipticity measurements cannot be attributed to a single com- 
ponent. For simplicity, we therefore reduce the difference of the 
ellipticity measurements to a one-dimensional quantity: 

Ae = -y/(ei,hst - £i,i) 2 + (fi2,hst - £2,i) 2 (B.l) 

Ellipticity measurements are considered to be equivalent if 
Ae < 0.2 and inconsistent if Ae > 0.2. In the first sample, 
there are 185 galaxies, whereas in the second there are 165. 

B.2.2. Magnitudes 

For the lensing analyses, the galaxy catalogs were split into 
three parts according to brightness to investigate their contri- 
bution to the lensing signal. The division was chosen such that 
each sample contained an equal number of galaxies. With the 



merged catalog, we can now examine how these samples com- 
pare in the two datasets. 

Fig. lB.3l illustrates the magnitudes measured of the matched 
objects. Most of the objects found in the /-band image are con- 
sidered "bright" objects in the HST data. The fact that we also 
detected tangential alignment in the bright HST bin confirms 
the lensing signal seen in the ground-based image. 

As would be expected, the shape measurements agree best 
for the brighter objects, while they are inconsistent for most of 
the faint objects. 

B.2.3. Pf correction 

One quality attribute we used previously to classify objects was 
the factor Pf by which the smearing due to the PSF was cor- 
rected. In Fig. IB.4l we plot an r g - Pf diagram of the objects 
found in the /-band image within the field covered by the HST 
mosaic. There is only a slight indication that ellipticity mea- 
surements with Pf < 0.4 are less reliable than others. 

B.2.4. Signal-to-noise ratio and weighting 

Two other criteria for the reliability of a shape measurement 
are the signal-to-noise ratio of a detection and the weight it was 
assigned. The latter one was taken to be inversely proportional 
to the variance <x 2 of the ellipticities of a galaxy ensemble with 
similar noise properties [see Sect. 12.3 .91 . Fig. IB. 51 shows the 
signal-to-noise ratios and weights of the /-band objects in the 
HST field. As expected, the ellipticity measurements deviate 
for those objects with a low signal-to-noise ratio. These are 
in general down-weighted, although the weight itself is not a 
clear indicator of the reliability. The weighting scheme could 
therefore be improved. 
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The lensing signal in the ground-based data was carried by 
faint galaxies with a signal-to-noise ratio of less than about 4. 
These are precisely those galaxies for which the space-based 
measurements give different ellipticities. It is therefore not sur- 
prising that the amplitude of the M ap signal is different in the 
two datasets. 

B.2.5. Dependence on chip position 

As mentioned before, we cannot judge the quality of the 
anisotropy correction of the HST images due to the small num- 
ber of stars in the images. With the correlated catalogs, we 
have another test of the correction: if it is faulty, then we ex- 
pect to see a systematic variation of the ellipticity differences 
Ag] = (e^hst - £1,1) and AE2 = (£2,hst _ £ 2,i) between the three 
chips and/or with position on the chip. We reduce the latter 
one to a one-dimensional quantity by considering the distance 
r from the respective chip center; this is motivated by the ob- 
servation that the anisotropy seen in the HST chips is largest at 
the edge (see Fig. lA.4> . 

In Fig. GO] we plot As 1 and As2 as a function of r, where 
we distinguish between the three chips. Clearly, an object's po- 
sition in the HST mosaic is not the main cause of discordance 
in the ellipticity measurement. However, there are trends visi- 
ble such as a slight overestimation of ei for small r in Chip 4 
and an underestimation of e\ for small r in Chip 3. To quantify 



3 



5 



Fig. B.l. The magnitudes and half-light radii of the /-band ob- 
jects positioned in the field observed by the HST. Open circles 
denote objects for which one HST counterpart was found, solid 
ones no counterpart, and crosses more than one counterpart. 
For some of the bright, large objects no counterpart was found 
due to the large masked areas in the HST mosaic. 
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Fig. B.3. Magnitude measured in the HST data vs. magnitude 
measured in the /-band image for the matched objects. The 
symbol coding is identical to that of Fig. IB. 21 The dashed lines 
indicate the cuts used to divide the respective catalogs into 
three samples. 
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Fig. B.6. Ellipticity differences hs\ and As2 as a function of position in the HST mosaic. Plotted on the x-axis is the distance 
from the chip center (measured in pixel units after the drizzle process); objects detected on Chip2 are shown as circles, Chip3 as 
triangles, and Chip4 as crosses. The distinction we make between equivalent and inconsistent shape measurements is illustrated 
by the black, long-dashed lines (if for the respective other component Ae, = 0). The dotted, short-dashed, and dot-dashed lines 
give linear fits to the respective data. 
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Fig. B.4. The ellipticity correction factor Pf and half-light radii Fig. B.5. The signal-to-noise ratios snr and weights of the I- 

of the /-band objects within the HST field. Crosses denote ob- band objects in the HST field. The color-coding is identical to 

jects for which no counterpart was found; the other colors are that of Fig. lB~4l 
the same as in Fig. IB. 21 



these, we fit a linear function Ae,- = m, • r/1000 + £>,, for each 
chip. We weight the Ae, values by the same weighting scheme 
as before, but do not assign error bars to the distance. The re- 



sults are: 



nt2 
b 2 



Chip 2 

-0.21 ±0.17 
0.14 + 0.10 



0.05 d 
-0.01 



0.17 
zO.10 



Chip 3 

0.23 + 0.16 
-0.13 + 0.10 

0.02 + 0.16 
-0.02 + 0.10 



Chip 4 

0.14 + 0.19 
-0.05 + 0.11 

0.03 + 0.19 
-0.03 + 0.11 
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For the S2 component, the ellipticity difference is remark- 
ably constant at zero for all chips. However, for the s\ com- 
ponent, there is a significant deviation from such a behavior in 
Chips 2 and 4. The anisotropy pattern of the HST chips is ap- 
proximately circular in appearance (see Fig. IA.4L so that an in- 
sufficient anisotropy correction should affect both components 
equally. The deviation in ei may therefore be a hint that the 
CTE and / or the drizzle process affect the ellipticities system- 
atically. On the other hand, the linear fits are certainly noisy 
due to the small number of objects, particularly at low r. 

Considering that Ae2 is on average constant at zero, but 
also that Aei is consistent with zero for a large number of ob- 
jects (400 < r < 800), strengthens our assumption that the 
anisotropy correction worked properly. 

B.3. Leasing analysis 

Since the ellipticity measurements agree on average, one can 
assume that the lensing analysis of our set of matched objects 
should yield similar results. However, the individual analyses 
indicate otherwise: the faint galaxies in the ground-based data, 
which caused most of the lensing signal, correspond to the 
bright (and medium bright) galaxies in the HST data, which 
gave only a weak signal. 

In order to directly evaluate the correlation (or discrepancy) 
between the ellipticity measurements and the lensing analysis, 
we perform several M ap analyses of the matched galaxies (see 
Table IB. ll for a quantitative summary). 

As a reference, we perform a M ap analysis of the 350 
matched objects using once the ground-based ellipticity mea- 
surements and once the space-based measurements (Fig. IB.7t . 
It is quite remarkable how well the M ap contours agree for 
these; however, for the dark clump peak, the space-based val- 
ues reach only half the height of the ground-based. 

Restricting the galaxy sample to those for which Ae < 0.2 
further highlights the agreement, as is expected (left hand side 
of Fig. lB~8l . These are mainly the galaxies considered "bright" 
also in the ground-based image (see Fig. lB.3> . which contribute 
to the lensing signal found there. Even for these, the M ap values 
of the dark clump peak are a factor of 1 .4 higher for the ground- 
based ellipticity measurements. 

For the 120" scale, the discrepancy in the M ap measure- 
ments at the dark clump position is even larger for those objects 
with Ae > 0.2 (right hand side of Fig. lB~8l . In the ground-based 
data, there is a > 3.7cr peak slightly to the right of the dark 
clump position, while the M ap values in the space-based data 
are comparable to the Ae < 0.2 sample, with a significance 
of ~ 2<x. At the dark clump position, the amplitudes in the 
ground-based values are on average a factor 3 larger than the 
space-based ones for all filter scales. Interestingly, for the 
smallest two filter scales, the highest significance found in the 
vicinity of the dark clump is equivalent in the two datasets. 
However, this is based only on a very small number of galaxies. 

As there is a M ap peak in all these samples, we can be sure 
that there is some tangential alignment around the dark clump 



candidate. But it is still unclear why it is measured to be so 
much larger in the ground-based image. 

B.3.1 . Interchanging the weights 

For the lensing analysis, the entries in the ground-based and 
space-based catalogs differ mainly in the shape measurements 
and in the weights assigned to each object. We have shown that 
the ellipticity measurements agree on average, but we have not 
yet considered the different weights assigned. 
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Fig. B.9. A comparison of the weights assigned to the matched 
objects in the HST catalog and the ground-based catalog. 
Again, solid points denote objects with Ae < 0.2, open circles 
Ae > 0.2. 



In Fig. IB. 91 we directly compare the weights assigned to the 
matched objects. Certainly, objects which were down-weighted 
in the ground-based data have higher weights in the HST cata- 
log. But else, the scatter is fairly large. 

To test the influence of the weights on the lensing analyses, 
we perform an analysis where we interchange the weights, i.e. 
we assign to each object in the ground-based image the weight 
of its counterpart in the HST image and vice versa. The results 
are listed in Table lB.ll It is interesting to note that at most filter 
scales, this causes the ground-based value to decrease and the 
space-based value to increase. But the M ap values are within the 
error bars of the ones with the original weights and are there- 
fore only slightly dependent on the weights. 

B.3.2. Recalibrating the HST ei component 

We noted earlier that the linear fit applied to the ei components 
of the matched objects yields ay-intercept with almost lcr sig- 
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Fig. B.7. M ap analysis of the 350 objects common to both the ground- and the space-based catalog. The upper two panels show 
M ap contours, starting at 0.02 and increasing in 0.02 increments, where a filter scale of 80" was used for the left panel, and 
a filter scale of 120" for the right one. The lower two panels give signal-to-noise contours, starting at 1.5cr and increasing in 
0.5cr increments. The white contours correspond to the analysis using ground-based ellipticities, and black ones to the one using 
space-based ellipticities. 



nificance (Sect. IB.2. II . This might point to a problem of the 
HST data related to its CTE or the drizzle procedure. 

Earlier, the ellipticity measurements from all three chips 
were considered. The ellipticities compared were those defined 
with respect to the rectascension axis. However, any problem 
regarding the CTE would bias the ellipticity measurements 
with respect to the read-out direction, which is different for 
all three chips. Also, the effect of the anisotropy correction is 
likely to be different for each chip. 

Therefore, we repeat the ellipticity comparison for the 
matched objects separated according to which HST chips they 
were measured in. The linear fits yield (with the same definition 
of the fit-parameters as in Sect. IB.2Tt : 



by 
b 2 



Chip 2 

1.09 ±0.31 
0.035 ± 0.047 

0.97 ± 0.23 
0.015 + 0.045 



Chip 3 

0.89 ±0.21 
-0.006 ± 0.043 

1.01 ±0.25 
-0.010 ±0.047 



Chip 4 

0.83 ± 0.22 
0.032 ± 0.047 

1.07 ±0.32 
-0.015 ± 0.054 



For the s 2 component, the agreement is still very good for all 
three chips. For the si component, the deviations from a slope 
of unity are larger. Particularly, Chip 2 and 4 exhibit a notable 
non-zero y-intercept. If these deviations are indeed a result of 
systematics in the HST data, applying the inverse relation to the 
HST measurements should then on average retrieve the ground- 
based ellipticities. 

We modify the HST measurements by the relation 
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for each chip i. We do not alter the s 2 components, as we deem 
their agreement with the ground-based data satisfactory. 

The results of a lensing analysis of the matched objects with 
these modified ellipticities are listed in Table IB. 11 Indeed, this 
transformation yields M ap values 25% larger than the original 
ones for the smallest two filter scales. For the 100" filter scale, 
the effect is smaller, and for the 120" filter scale it is negligible. 

If the height of the M ap peak were uncorrelated with the 
systematic deviation in the ellipticity measurement, the modi- 
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Fig. B.8. M ap analysis of objects with As < 0.2 (left side) and Ae > 0.2 (right side). The upper panels show the M ap values 
measured in the 120" filter scale, the lower ones the corresponding signal-to-noise ratios. The color coding and the contours are 
identical to Fig. IB. 71 



fication would only add noise to the M ap statistic, as most ellip- 
ticities are amplified due to m, < 1 . The M ap value of the peak 
therefore would not be significantly altered, as is the case at 
least for the largest filter scale. For the other filter scales, there 
is a slight increase in both M ap and SNR value. But the change 
is at most 0.3cr, so the variations are still within the standard 
deviation of the original measurement. 

The results of a lensing analysis of the complete HST cat- 
alog with this modification of the ellipticity are very similar to 
the original one. But as the linear fit which was the basis of this 
modification applies only to bright objects in the HST catalog, 
it is very speculative to extrapolate this to fainter objects. 

As this modification alone is not able to reproduce the large 
lensing signal seen in the ground-based data, particularly at the 
120" filter scale, the offset of the space -based ellipticities is not 
the cause of the discrepancy in the M ap measurements. 

B.3.3. Rebinning the galaxies in the HST catalog 

For the lensing analysis of the HST data, we had divided the 
galaxy sample into three magnitude bins of equal numbers of 



object. About half of the objects in the brightest bin were also 
detected in the /-band image and are used in the comparisons 
presented in this chapter. For these, we have confirmed the 
presence of tangential alignment, even though the M ap ampli- 
tude differs in the two datasets. But we fail to detect alignment 
in the medium bright HST bin. 

The completion limit of the ground-based image falls 
within the brightest HST bin (cf. Fig. IB.3i . and so a number 
of objects in the latter, though of similar brightness, were not 
detected in the ground-based image. To test how much the ad- 
ditional objects in the HST image contribute to the lensing sig- 
nal, we rearrange the HST brightness bin: instead of the mag- 
nitude cut between the brightest and medium bin, we split the 
galaxies according to whether or not they are a counterpart of 
an object detected in the /-band image. From the latter sample 
the faintest bin is split with the same magnitude cut as before. 
This effectively moves several objects from the bright into the 
medium bin and only a few the other way. The faint bin remains 
almost unaltered. 



In Fig. IB. Ill we present the signal-to-noise contours of a 
M ap analysis of these three bins. The medium bin, i.e. galaxies 
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Table B.l. Overview of the M ap analyses discussed in this appendix. Shown are the results at the position of the dark clump 
as measured in the ground-based data, since we want to determine the cause of the high M ap amplitude in the ground-based 
data at this position. We list the ellipticity dispersion, the M ap value itself, and the signal-to-noise ratio. For those analyses of 
the 350 matched objects, we give both the results of using the ground-based and of using the space-based data. For the second 
block from the bottom of the table, the ground-based analysis uses the space-based weights and vice versa. For the last block, the 
si component of the HST measurements is modified by applying the inverse of the linear fit between space-based and ground- 
based ellipticities. Evaluating M ap only at the ground-based dark clump position biases the result in favor of larger ground-based 
values. For a better comparison of the general values, we quote the highest significance found in the vicinity of this position in 
parentheses. 
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brighter than m = 26.5 that were not detected in the ground- 
based image contains the most galaxies (835 compared to 350 
in the bright bin and 594 in the faint bin), so that the noise is 
smallest. Yet, the M ap values in the vicinity of the dark clump 
are compatible with zero within lcr. 

Apparently, there is no tangential alignment present in the 
objects that were moved from the bright bin to the medium bin. 
This is compliant with the observation that for those objects 
also detected in the /-band image, the signal-to-noise ratio of 
the M ap peak is larger than for the bright magnitude bin (cf. 

Fig. CD- 

Intriguingly, this implies that of the faint galaxies at the 
completion limit of the /-band image, those with a tangen- 
tial alignment to the dark clump were preferentially detected. 
However, the fact that there is no tangential alignment for about 
half of the objects in the HST image, which are neither at the 



bright nor the faint end of the magnitude distribution, clearly 
disfavors the hypothesis that the tangential alignment in the 
brighter galaxies is caused by a lensing mass. 



B.4. Comparison of the tangential ellipticities 

So far, we have not found any systematic that could be the sin- 
gle cause for the discrepancy in the lensing analyses. However, 
it is clear that for the faint galaxies which mainly caused the 
lensing signal in the ground-based image, the ellipticities agree 
with the HST measurements only on average, but not on an 
object-to-object basis. This is not surprising as the noise in the 
ground-based image is likely to influence the shape measure- 
ments. Yet, it is unlikely that noise can cause (or amplify) tan- 
gential alignment around a certain point. 
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Fig. B. 10. Comparison of the ellipticity measurements in the HST image and the /-band image, separated by HST chip. The 
upper row shows the s\ measurements, the lower the S2 measurements. Filled circles denote objects with As < 0.2, open ones 
As > 0.2. A weighted linear fit is indicated by the solid line. 
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Fig. B.12. Comparison of the tangential ellipticities measured 
toward the dark clump (as measured in the ground-based data) 
of space-based and ground-based objects. Objects with Ae < 
0.2 are shown as solid, objects with Ae > 0.2 as open circles. 
The dashed line indicates a linear fit to the first sample, the 
dotted line one to the latter sample. The solid line represents 
the fit to the combined set. 



In Fig. lB.12l we compare the tangential ellipticity s t of the 
matched objects with respect to the dark clump centroid found 
in the ground-based data. Only objects within 120" of that po- 
sition are shown. Evidently, most objects in the ground-based 
data have a positive e t . The scatter is comparable to the com- 
parison of the individual ellipticity components (Sect. IB.2.U . 
Likewise, we fit a linear function, where we distinguish be- 
tween the complete sample of matched galaxies and how well 
the ellipticity measurements agree (as before, the fit is ex- 
pressed such that the ground-based data are the independent 
variable): 

all Ae < 0.2 Ae > 0.2 

m 0.98 ±0.21 1.03 + 0.28 0.93 + 0.32 

b -0.047 + 0.043 -0.018 + 0.055 -0.086 + 0.071 

For each of these fits, the y-intercept is negative, which 
indicates larger ground-based values. While it lies well within 
the error bars of the fit for those galaxies with Ae < 0.2, it 
is significantly non-zero for those with Ae > 0.2. For those, 
also the slope indicates a preferentially larger e t value in the 
ground-based data. 

On the basis of the correlation of the individual measure- 
ments of the tangential ellipticity, the agreement between the 
two datasets is comparable to that of the ellipticity components, 
as would be expected. Yet, the mean tangential ellipticity is dif- 
ferent: for the ground-based data, it is (e t j) = 0.086 + 0.017, 
for the space-based data, (e t ,h s t) = 0.034 + 0.017. These are 
very similar to the results of the M ap analysis at the same po- 
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Fig. B.ll. M ap analysis of the HST data, with the galaxies split into three samples according to whether they are a counterpart of 
an object detected in the /-band image (white contours). Those objects not detected in the /-band image are split by a magnitude 
cut at m = 26.5, where gray contours corresponds to galaxies brighter than this limit, black contours to those fainter. The contours 
start at 1.5<x and increase in 0.5cr increments. The filter scales are 80", 100", 120", and 140" (from upper left to lower right). 



sition at the 120" filter scale, M^ p (120") = 0.087 + 0.019 
and A/j^(120") = 0.043 ± 0.019 (compare Table ED- This 
again demonstrates that the coherent alignment is present over 
a range of distance from the dark clump, as the M ap statistics 
is a weighted average tangential ellipticity, where the weight is 
dependent on the distance. 

Within the standard deviations of the cited values, the two 
analyses are compatible with each other. However, these errors 
are not independent, as they are based on the ellipticity dis- 
persion of the same galaxies. To estimate the significance with 
which the ground-based mean is larger than the space-based 
mean, we perform a bootstrap analysis, i.e. from the 195 galax- 
ies in the sample, we draw at random 195 objects, with replace- 
ment. From these we calculate (e t ,i) and (e t ,hst) and repeat the 
procedure 1,000,000 times. The result is shown in Fig. IB.131 
We find that the space-based mean exceeds the ground-based 
mean in only 176 cases, i.e. the latter one is larger with 3.8<r 
confidence. 

This analysis is biased in the way that the reference posi- 
tion is the dark clump position as measured in the ground-based 
data. Therefore, we repeat the analysis for the centroid found 



in the space-based data. Even for this case, the ground-based 
mean value exceeds the space-based mean with 3.1cr confi- 
dence. 

It is quite puzzling that despite the general agreement of 
the ellipticity measurements, the mean tangential alignment 
can differ so much. As yet, we have no explanation for the 
cause of this. 

In Fig. IB.14l we present clippings of some of these galax- 
ies from the two images to convey a visual impression on how 
they compare in the images and how the shape measurements 
relate to the image, especially for the ground-based image. We 
find that in the HST image, the measured ellipticity very well 
represents the shape of the object. This is to be expected as 
all these objects are detected with a high signal-to-noise ratio 
in the HST image. In the /-band image, these objects have a 
low signal-to-noise ratio and the ellipticity measurements only 
vaguely reflect the shape of the objects. About 10% of the ob- 
jects are actually two or more objects which were not resolved 
in the ground-base image; for about 20% a nearby (though re- 
solved) object seems to influence the ellipticity measurement. 
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Fig. B.13. Results of a bootstrap analysis of the mean tangen- 
tial ellipticity of the matched objects within 120" of and with 
respect to the dark clump position in the ground-based data. 
For objects below the long-dashed line, (e t ,i) > (£t,hst)> above 
it, (s t j) < (e t ,hst)- The solid black (red) line and dashed black 
(red) line indicate (e, j) ((s t> h s t)) without bootstrapping. 



It is especially striking that for many, the ellipticity modulus |e| 
is clearly overestimated. 

B.5. Summary 

We have shown that on average, ground-based and space-based 
ellipticity measurements agree very well. This justifies the as- 
sumption that weak lensing analyses based on ground-based 
data yield reliable results despite the smearing of object shapes 
by the Earth's atmosphere. 

For the HST ellipticities, we have found only little evidence 
that they might be biased due to the CTE problem, the PSF 
unders amp ling, or the anisotropy correction. This holds at least 
for bright objects, for fainter ones we have no possibility of 
comparison. Yet, previous analyses (Sect. IA.2l have shown that 
these problems have little influence. 

The fainter an object, the less reliable its shape measure- 
ment is (Fig. IB.3t . The strong lensing signal seen in the ground- 
based data is caused mainly by faint objects, for which the 
ellipticity measurements in the two datasets deviate. In both 
datasets, there is a M ap peak at the position of the dark clump, 
but the degree of tangential alignment is much stronger in 
the ground-based data. For many of these objects, the elliptic- 
ity modulus |s| is overestimated in the ground-based measure- 
ments, and many of these seem to be accidentally tangentially 
aligned to the dark clump, thus causing a high M ap signal. 
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Fig. B. 14. Direct comparison of objects with Ae > 0.2 within a 120" radius around the dark clump (position measured in the 
ground-based data) in the /-band (left side) and HST image (right side). For each object, we list its ellipticity components ei and 
£2 in both images, which are also indicated by ellipses superposed on the images. The signal-to-noise ratio snr of each detection 
is listed as well. In the ground-based image, we also indicate the direction to the dark clump (by an arrow in the image) and its 
distance dist, and list the tangential ellipticity e t towards the dark clump for both images. For the HST images, we also give each 
object's distance d(i) to the center of the chip ;' it was measured on. Only 8 exemplary objects are shown. 
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